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Measurements of Baryonic Acoustic Oscillations in galaxy surveys have been recognized as a pow- 
erful tool for constraining dark energy. However, this method relies on the knowledge of the size 
of the acoustic horizon at recombination derived from Cosmic Microwave Background Anisotropy 
measurements. This estimate is typically derived assuming a standard recombination scheme; ad- 
ditional radiation sources can delay recombination altering the cosmic ionization history and the 
cosmological inferences drawn from CMB and BAO data. In this paper we quantify the effect of 
delayed recombination on the determination of dark energy parameters from future BAO surveys 
such as BOSS and WFMOS. We find the impact to be small but still not negligible. In particu- 
lar, if recombination is non-standard (to a level still allowed by CMB data), but this is ignored, 
future surveys may incorrectly suggest the presence of a redshift dependent dark energy compo- 
nent. On the other hand, in the case of delayed recombination, adding to the analysis one extra 
parameter describing deviations from standard recombination, does not significantly degrade the 
error-bars on dark energy parameters and yields unbiased estimates. This is due to the CMB-BAO 
complementarity. 



I. INTRODUCTION 

A key goal of modern cosmology is to investigate the 
nature of the dark energy component, responsible for 
the current accelerated expansion of the Universe. A 
variety of obscrvablcs, from Cosmic Microwave Back- 
ground (CMB) anisotropies [2, Q to galaxy surveys 
H> BL IE IE 13 1 wm continue to be measured with in- 
creasingly refined accuracy in the next years thanks 
to a coordinated effort of satellite, ground based and 
balloon-borne missions. Despite the fact that it ac- 
counts for about 70% of the total energy density of 
the universe, dark energy is largely unclustered and is 
typically measured just by its effect on the evolution 
of the expansion history (i.e. the Hubble parameter) 
at redshift z < 2. Since the cosmic expansion de- 
pends on other key parameters as curvature or matter 
density, the nature of dark energy can therefore be 
revealed only by combination of different observables 
and/or observations over a wide redshift range. 

A recent development in the field of large-scale 
galaxy surveys, is to use the measurement of Bary- 
onic Acoustic Oscillations (BAO) [1, [E [l(| signal at 
z < 2 as a standard ruler. Oscillations in the pri- 
mordial photon-baryon plasma leave an imprint in the 
matter distribution. Since the frequency of these os- 
cillations is related to the size of the sound horizon 
at recombination, which is well constrained by CMB 
measurements, it is possible to use the measurements 
of those oscillations at different redshifts as a standard 
ruler and therefore determine the rate of the cosmic 
expansion. 

It has been found that a number of theoretical sys- 
tematics including non- linear growth, non- linear bias, 
and non-linear redshift distortions can be efficiently 



modeled to minimize their contribution to uncertain- 
ties in the analysis of BAO observations [ll|, EH, EH ■ 
As such, they are one of the key observables of the 
next decade in cosmology and it is therefore impor- 
tant to investigate how solid are the theoretical as- 
sumptions behind this method. A crucial assumption, 
indeed, is the possibility of having a very accurate 
meausurement of the size of the acoustic horizon at re- 
combination. fl7T ] show that the interpretation of low 
redshift BAOs is robust if CMB correctly determines 
the baryon-to-photon ratio and the matter-radiation 
equality but not the matter baryon or photon densi- 
ties alone. CMB measurements of these quantities is 
generally robust, but , in view of high-precision data, 
it is important to quantify any possible systematic ef- 
fects introduced by deviations from the standard evo- 
lution of the early universe . Previous papers have 
studied, for example, the stability of the result under 
the hypothesis of an extra background of relativistic 
particles or early dark energy ([131). In this case, the 
epoch of matter-photons equality is shifted and this 
may bias the determination of size the acoustic hori- 
zon if unaccounted for. 

In this paper we consider another possible mecha- 
nism that could, in principle, affect the BAO inter- 
pretation as standard rulers if not modeled, namely a 
possible delay, or change in the recombination epoch 
that would result in a variation of the acoustic horizon. 
Several papers, rece ntly , have indeed considered this 
hypothesis (see e.g. [3 IHI, EE [H, H3, [HI, H3] ) - Dark 
matter decay or annihilation, black hole evaporation, 
or cosmic string decay can, for example, produce an 
extra background of resonance photons that could de- 
lay recombination (see and references therein). 
More exotic mechanisms as, for example, variations 
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in the fine structure constant or in the Newton con- 
stant, can also lead to a modified recombination epoch 
([Ullill])- We conclude that current constraints on de- 
layed recombination are compatible with a change in 
the sound horizon which is non-negligible when com- 
pared with the precision for future BAOs data. This 
means that, if unaccounted for, delayed recombination 
could bias dark energy constraints from BAO. On the 
other hand, if delayed recombination is accounted for 
in the joint CMB-BAO analysis form forthcoming sur- 
veys, error-bars on dark energy parameters are only 
degraded by less than 10 % compared to the standard 
analysis. 

The paper is organized as follows: in the next sec- 
tion we discuss delayed recombination and its effect on 
the sound horizon. In Sec. Ill we review the Fisher 
matrix formalim, while in Sec. IV we forecast the 
impact of delayed recombination on dark energy pa- 
rameters inference for several future surveys. Finally, 
in Sec. IV, we draw together our findings and discuss 
the conclusions and their implications for future work. 



where e a is assumed constant with redshift. 

This leads to a modified evolution of the ionization 
fraction: 



where std refers to the the standard scenario. 

The introduction of e a has the main effect of delay- 
ing the redshift of recombination z*. In the standard 
scenario this redshift is given by ([27j]): 

z* = 1048(1 + 0.00124^°- 738 )(l+gi< 2 ) (5) 
,9i = 0.0783Lu h °- 238 /(l + 39.5w£- 763 ) 
.92 = 0.560/(1 + 21. l^- 81 ) 

where io m and uib are the matter and baryon phys- 
ical energy densities respectively. Including positive 
values of e a decreases z* as (PJl): 



II. DELAYED RECOMBINATION AND THE 
SOUND HORIZON 

In this section we modify the standard recombi- 
nation picture by including the possibility of extra 
sources of resonance (Lyman-a) photons. Note that 
we do not modify matter-radiation equality, an effect 
already explored by [l7| ■ 

In the standard recombination scenario (see [25l . 
l2l| ) the evolution of the electron ionization fraction, 
x e is given by: 



z*(e a ) = {l + 3e a )~" MA z,(e a = 0) 



(6) 
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where x e is the electron fraction per hydrogen nuclei, 
77 is the number density of atoms, a c and b c are the 
effective recombination and photo-ionization rates for 
principle quantum numbers > 2, AB is the difference 
in binding energy between the I s * and 2 nd energy lev- 
els and 



C = 



1 + A'A ls2s 77 ls 
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K = 



8nH(z) 
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where X a is the wavelength of the single Ly-a transi- 
tion from the 2p level, Ai S 2 S is the decay rate of the 
metastable 2s level, n\ s = n(l — x e ) is the number of 
neutral ground state H atoms, and H(z) is the Hubble 
expansion factor at a redshift z. 

As in [3, EH, [l(| , we include the possibility of ex- 
tra resonance (Ly-a) photons at recombination with 
number density n a which promote electrons to the 2p 
level (see 0, M, M EH): 



dn a . . 

— — = e a H{z)n 
at 



(3) 



Conservative limits on e a from current CMB data 
are of the order of e Q < 0.5 at 95% c.l. The 
redshift of recombination can therefore be smaller by 
^5% respect to the standard case. These constraints 
were obtained for a LCDM cosmology: the bound on 
e a would be relaxed even further for models where 
dark energy is not a cosmological constant. 

We can easily understand the main effect of a de- 
layed recombination by looking at the variations in the 
visibility function which describes the density proba- 
bility for last scattering for a photon at redshift z. Its 
peak defines therefore the epoch of recombination. As 
we can see from Figure [H where we plot the visibility 
function defined as in [29( for different values of e a , 
an increase in e Q shifts the visibility function towards 
lower rcdshifts and increases the width of the distribu- 
tion. We can therefore expect two main effects on the 
CMB anisotropics: a displacement of the peaks in the 
anisotropy and polarization power spectra due to the 
shift in the recombination epoch and a simultaneous 
damping of the anisotropies (see Fig. [2]) due to the 
increase in the finite thickness of the last scattering 
surface. 

In what follows wc quantify the impact of this effect 
on the acoustic oscillations imprinted in the matter 
power spectrum that can then be used as a "standard 
ruler" at lower redshifts. We investigate the variation 
of the sound horizon size s as a function of z* and, 
therefore, e a . In fact galaxy surveys, by measuring 
the scale along and across the line of sight, constrain 
H(z)~ 1 /s and Da(z)/s, where Da is the angular di- 
ameter distance at redshift z. It is therefore straight- 
forward to expect that a systematic change in s could 
bias the determination of iJ(z) -1 , Da{z) and the de- 
rived parameters. 
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The size of the sound horizon s at recombination 
can be approximated by the following formula ([2S|]): 



V3H 
19.8 Mpc 



da 



hi 



'o y/(a + a cq ))(l + R{a)) 

' V R(a*) + R(a cq ) + y/l + R{aS 



1 + v/i?(a cq ) 



where Hq is the Hubble constant, a is the scale fac- 
tor normalized to unity today, R(a) = 30496a->bi, the 
equality is at a cq = l/(24185w m ) and a* = + 

Shifting the epoch of recombination, while leav- 
ing as unchanged the remaining parameters, clearly 
changes the extrema of the integral and the size of 
the sound horizon. 

In Figure [3] we show the percentage variation in ,s 
as a function of e a : a value of e a as large as 0.5 yields 
a few percent variation in the CMB estimate of the 
sound horizon. This can potentially propagate into a 
large shift in the recovered dark energy parameters: 
at z ~ 0.6 a ~ 1% error in s yields a ~ 4% error 
in w assumed constant. This shift is comparable to 
the statistical error on this quantity expected from 
forthcoming surveys. However, as we will see in the 
next sections and already from Figure [21 the impact of 
delayed recombination on obscrvablcs is broader. In- 
deed, the damping in the anisotropy power spectrum 
introduces a degeneracy with other parameters as the 
scalar spectral index n s . As showed in [l6| consider- 
ing values of e a ~ 0.1 could mimic variations as large 
as ~ 10% in n s and in several other variables as the 
Hubble constant. 



III. FISHER MATRIX ANALYSIS FOR BAO 
AND CMB: STANDARD RECOMBINATION 

In this section we briefly review the Fisher formal- 
ism mostly based on the seminal paper by [30l | and we 



refer the reader to this work for further clarifications. 
At the end of this section we cross-check our results 
with similar analyses already presented in the litera- 
ture and obtained under the assumption of standard 
recombination. 

The usual definition of the Fisher matrix is (see 
[Mil): 



Fa/3 = (- 



d 2 InL , 
dp a dpp ' 



(8) 



where L is the likelihood function for a set of param- 
eters pi. If we suppose that the likelihood has a max- 
imum in some point of the parameter space p®, the 
fiducial model, and if the data have a normal distri- 
bution then the inverse of the Fisher matrix represents 
the covariance matrix of the parameters, i.e.: 



((p a -p° a )(p ! 3-p o )) = (F- 1 U, 



(9) 



as a consequence the statistical errors on the param- 
eters are given by the square root of the diagonal ele- 
ments of the inverse Fisher matrix. According to the 
Cramer- Rao inequality ([IH) for unbiased estimators 
this is the best statistical error that one can obtain 
for the generic parameter p a : 



<J Pa > ^{F~ l ) aa . 



(10) 



For a galaxy survey the Fisher matrix can be ap- 
proximated as (Hi||): 
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FIG. 1: Visibility function (defined as the density prob- 
ability of photon last scattering), in function of e a . The 
inclusion of a delayed recombination shifts the peak of the 
distribution at lower redshift and increases its width. 



3500 

3000 

=1 2500 



q — 2000 

IF 

CM 1500 



_+ 1000 
500 











— ta - ■ 










— ea - 0.1 










ea - 1 










ea - 10 


1 


\ '-■/ 


\ s \ V'* x ^ \ 
1 1 


i 




200 


400 


600 800 


1000 


1200 1400 



FIG. 2: Angular power spectrum of CMB anisotropies in 
function of e Q . The inclusion of a delayed recombination 
shifts the peaks of the spectrum and damps the structure 
at small angular scales. 
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where the derivatives of the linear matter power spec- 
trum P{k 1 ^x) are evaluated at p° and V e ff is the ef- 
fective volume of the survey, given by: 



V eff {k,n) 



n g {f)P{k,pt) 



_n g (r)P(k,n) + 1 
n g P(k,fj) 



dr 



n g P(k,fi) + 1 



(12) 



V„ 



where in the last equality is assumed that the comov- 
ing number density of galaxies n g is constant. The 
quantity /i is the cosine of the angle between the unit 
vector along the line of sight f and the wave vector 
k, n = k • f/k. The integration over k that appears 
in (jlip is performed only up to a k max , to exclude 
non-linear scales. The value of k max is redshift de- 
pendent and we calculate it using the same criterion 
of [3(|, i.e. requiring a(R) = 0.5, where R = n/2k. 
The maximum scale of the survey, corresponding to 
kmim doesn't influence the results and we can assume 

^min 0. 

We consider two future surveys: the Baryon Os- 
cillation Spectroscopic Survey (BOSS [3f|), that will 
cover redshifts up to z = 0.7, and the Wide-Field 
Multi-Object Spectrograph (WFMOS (37|) splitted 
in two samples, the first (WFMOS1) at redshifts 
0.5 < z < 1.3 and the second (WFMOS2) covering 
2.3 < z < 3.3. We therefore divide the redshift space 
in six bins, one up to z = 0.7 for BOSS, four bins 
for WFMOS1 and one bin for WFMOS2. In Table 
U we list the survey area and the number of galaxies 
assumed in the analysis. 

Following [30l | , we consider 6 cosmological parame- 
ters, the physical matter density fl m h 2 , the physical 
baryon density flbh 2 , the matter fraction f2 m , the op- 
tical depth to reionization r, the scalar spectral index 
n s and the normalization A s . The fiducial model is 
given by Q m h 2 = 0.1274, Q b h 2 = 0.021, Q m = 0.26, 
t = 0.05, n s = 0.96, A s = 2.4 ■ 10" 9 and h = 0.7. We 




FIG. 3: Percentage of variation in the size of the acoustic 
horizon at recombination in function of e a . Delayed re- 
combination with e a ~ 0.5 could increase the horizon by 
few percent. 
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TABLE I: Experimental specifications for the surveys used 
in this paper showing redshifts, survey area (A) and num- 
ber of galaxies observed (N). 



always assume flatness, i.e. £7 a = 1 — I n addi- 
tion to these parameters we included other, redshift- 
dependent, parameters as the angular diameter dis- 
tance Da(z), the Hubble parameter H(z), the linear 
growth factor G(z) and the linear redshift distortion 
j3{z). Our fiducial values for H(z) and Da (z) are ob- 
tained assuming a cosmological constant model. For 
the bias we assume b = 2 at z = 0.35, b = 3 for WF- 
MOS2 at z = 2.8 and values ranging from b — 1.2 at 
z = 0.5 to b = 1.7 for z = 1.3 for WFMOS1 With this 
parameterization each bin in redshift has a total of 10 
parameters, 6 common to all bins and 4 different for 
each bin. 

The observed power spectrum is given by: 



(13) 

where the first factor accounts for the fact that the 
observed cosmology (subscript r) could differ from 
the "true" cosmology, b is the bias, given by b — 
^m(z) ' 6 / P(z), and the term (1 + /3/x 2 ) 2 describes the 
linear redshift distortions. The linear growth factor 
G(z) is given by the ratio d(z)/6(Q) between the lin- 
ear density contrast at redshift z and at z = 0. We 
calculate the linear matter power spectrum at redshift 
z = 0, P°(z = 0, k), using the numerical code CAMB 
[H| . For convenience we take the reference cosmology 
to be equal to our fiducial model. 
The resulting Fisher matrix is combined with the 
CMB information from the Planck satellite [3^]. For 
a CMB experiment the Fisher matrix is given by [36j : 



rpCMB 
a/3 



1=2 i,j 



(14) 



where the C\ 3 are the well known power spectra for 
the temperature (TT), temperature-polarization (TE) 
and E-mode polarization (EE) (i and j run over TT, 
EE,TE) and Cov; is the spectra covariance matrix. 

We use the experimental configuration as described 
in Table HT1 with Z max = 1500 to calculate the sum in 
(fT4"|) . For the CMB Fisher matrix we include as a free 
parameter also the angular diameter distance to the 
last scattering surface Da,cmb- 

The full Fisher matrix is then simply given by the 
sum of the BAO Fisher matrices at each redshift bin 
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TABLE II: Experimental specifications for the Planck 
satellite. Channel frequency is given in GHz, FWHM in 
arcminutes, and noise in 10 -6 . 



Parameter a 
Q. m h 2 0.0007 
fi A 0.013 
w 0.11 
wi 0.14 

TABLE III: Uncertainties on dark energy and cosmological 
parameters from our Fisher matrix calculation in the case 
of standard recombination. 

and the CMB: 

N 

F t„ t = F CMB + S^ F! . (15) 

z=l 

For TV bins the total number of parameters is then 
6 + N x 4 + Ncmb where Ncmb is the number of addi- 
tional parameters for the CMB, in this case Da,cmb 
(and hence Ncmb = 1)- In the standard recombina- 
tion case we have 6 bins in rcdshift and our full Fisher 
matrix has 31 parameters. 

As explained in [3(j, by taking Da(z), H(z) and G(z) 
as separate parameters at each redshift one is avoiding 
any assumption about a specific dark energy model. 
The main advantage of this procedure is in consider- 
ing a Fisher matrix which is independent on the dark 
energy model and that can specified afterwards with 
a simple change of parameters. Here we allow for a 
possible time-dependent equation of state assuming a 
linear evolution with redshift: 

w(z) = wo + w\z (16) 

To obtain constraints on dark energy we first select 
the set of parameters that affect distances and depend 
on dark energy, i.e. p = (fl m h 2 ,Q m , Da{z)' s, H(z)' s) 
and we then marginalize over the remaining "nui- 
sance" parameters. We implement this by inverting 
Ffj* and selecting from the inverted matrix the ele- 
ments that correspond to the parameters of interest. 
We then invert this submatrix to obtain a Fisher ma- 
trix (Fj lm ) just for fl m h 2 ,Q m , D A {z)'s and H(z)'s. 
Finally we project in a new parameter space, for ex- 
ample q = (n m h 2 ,flA,wo,wi), by equating the log 
likelihood functions in the old and in the new param- 
eter space: lnL(p) = lnL(q). In the Fisher matrix 



formalism this corresponds to a contraction of the old 
reduced Fisher matrix F' nm with derivatives between 
old parameters and new parameters: 

_ dPn „, dp 
ij ~ ~dq~ nm ~dq~ [ 1 

The fiducial model is specified in the calculation 
of partial derivatives that appear in (|17[) . we there- 
fore assume a A-CDM model given by w = —1 and 
u>i =0. The results from our Fisher matrix analy- 
sis for the dark energy parameters are shown in Ta- 
ble IIIII The results, obtained under the assumption 
of standard recombination, are substantially in agree- 
ment with those reported in (30j . with errors on wo 
and w\ of a factor ~ 2 smaller since we are consid- 
ering improved experiments, with higher survey area 
and volumes. An improved Fisher matrix formalism 
has been recently presented in [j^ „ based on numer- 
ical simulations of [ifjj and where nonlinear effects, 
bias and redshift distortions are more accurately im- 
plemented. The errors on Da and H from their pub- 
licly available code for the BOSS survey at redshift 
z = 0.35 are respectively ~ 1.1% and ~ 2.1% to be 
compared with ours ctd/D ~ 1.5% and ajj /H ~ 2%. 
We can therefore conclude that, for the scope of our 
paper, the results are in reasonable agreement. 

Let us now describe in more detail the correlation 
between the different parameters in order to empha- 
tize the importance of the CMB measurement of the 
sound horizon. In the BAO only Fisher matrix, strong 
correlations are present between the angular diameter 
distances Da(z)'s and H(zys due to the poor deter- 
mination of £l m h 2 , flbh 2 and n s . These degeneracies 
can be reduced by an accurate CMB measurement of 
these cosmological parameters. As we can see from 
Table IIVI the correlations between the Da(z)'s and 
H{z)'s (calculated as p = {F- 1 ) lj /^/{F- 1 ) u (F- 1 ) jj ) 
from galaxy surveys alone arc large, showing that, 
until the scale of sound horizon and other parame- 
ters are unknown, only the product H(z)Da(z) can 
be precisely determined from BAO surveys. Is clear 
from the above discussion that if the CMB informa- 
tion is biased by the assumption of a incorrect model 
of recombination, it will have an impact on the de- 
termination of the cosmic parameters from BAO. We 
investigate this in detail in the next section. 



IV. DELAYED RECOMBINATION AND 
COSMOLOGICAL INFERENCES DRAWN 
FROM CMB AND BAO DATA 

As shown in [4l[ it is easy to compute the shift in the 
best fit parameters when other parameters are fixed 
at a wrong fiducial value within the Fisher matrix for- 
malism. In general if we fix a number p of parameters 
(say tpj, with 7 = 1, ...,p) to an incorrect value which 
differs from the true value for an amount Stp^, then 
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0.35 


8.39 


8.48 


-0.964 


1.48 


1.95 


-0.151 


0.6 


8.57 


8.82 


-0.890 


2.76 


3.48 


-0.160 


0.8 


8.16 


8.27 


-0.942 


1.91 


2.38 


-0.157 


1.0 


8.10 


8.17 


-0.957 


1.70 


2.07 


-0.212 


1.2 


8.01 


8.06 


-0.964 


1.57 


1.89 


-0.225 


2.8 


7.93 


7.24 


-0.979 


1.16 


1.29 


-0.217 



TABLE IV: Relative uncertainties on Da{z)'s and H(z)'s and their correlation coefficient (p) at each redshift. The 
correlation is large unless considering CMB observations, needed to calibrate the sound horizon. 



the others n parameters (0 a , with a = 1, n) will be 
shifted of an amount given by: 

60 a = -(F , - 1 ) af iSf h 6fy a,p = l,..,n (18) 

where (F'~ l ) a j3 is the inverse Fisher matrix in the 
space of 6 a parameters (an n x n matrix) and Sp^ is 
a subset of the full (n + p) x (n + p) Fisher matrix in- 
cluding also the tp a parameters. Equation (|18| shows 
that the shift on 9 a parameters depends on the corre- 
lations among various parameters, on how well these 
parameters are constrained and on how strong is the 
dependence of the observable from the p parameters 
fixed, i.e. on the derivatives of the observable with 
respect to the ipj parameters, which are contained in 
the Sp-y matrix. In our case we want to see what is the 
effect of assuming standard recombination, i.e. fixing 
e a = 0, while its true value is different from zero and 
we take it to be e a = 0.05. 

To study the effects of a modified recombination on 
BAOs survey we follow the method described above 
computing first a Fisher matrix for the parameters 
listed in the previous section and then repeating the 
computation including also e a . We then select from 
this second Fisher matrix the sub-matrix S^ 7 that ap- 
pear in (fl~5|) , i.e. the row and the column correspond- 
ing to e a , and calculate the shift for the distances, 
Hubble parameters and dark energy parameters. 

z SD/D SH/H 

0.35 -3.0% 3.0% 

0.6 -3.4% 3.3% 

0.8 -3.6% 3.6% 

1.0 -3.1% 3.0% 

1.2 -3.4% 3.3% 

2.8 -4.5% 4.5% 

TABLE V: relative shift in the measure of Da(z) and H(z) 
for the mean redhsift of each bin as computed from (|18|l 
assuming e a = in the hypothesis that true value is e a — 
0.05. 



We expect a shift 5s in the value of sound horizon to 
cause a shift in angular diameter distances and Hubble 
parameters, being Da oc s and H oc 1/s. Moreover is 
easy to verify that relative shift in Da(z) and H(z) 
at each redshift should be of the same extent and op- 
posite: where SDa/Da = —SH/H. This is shown 
in Table V. Applying (fl"5)) to calculate shift in Da(z) 
and H (z) confirms this, as shown in Table [V] We 
show also absolute shifts in Da and H for each bin 
in Figure SI The shifts in Da and H are at the 3%- 
4% level and, except for the last bin, are below the 
95%C.L. 

We remind that here we choose a fiducial model for the 
dark energy parameters with wq = — 1 and Wi = 0. 
Figure O shows that assuming an incorrect recombi- 
nation model introduces a shift in the recovered value 
of wq towards wq > — 1 and of w\ towards negative 
values, as a consequence of the shifts on Da and H . 
For Ae a = 0.05 the effect is comparable to the lcr 
error on wq, and slightly larger for w\. With the in- 
correct recombination model, the best fit are shifted 
to wo = —0.88 and w% = —0.22 respect to the fiducial 
A-CDM model. An incorrect assumption on recom- 
bination could therefore mislead us towards claiming 
deviations of dark energy from a cosmological con- 
stant. 

Even if this misinterpretation of the early universe 
physics induces shifts < 95%C.L. both for wq and 
u>i, these results show that an incorrect calibration 
of sound horizon due to a wrong assumption for the 
recombination model could bias estimation of dark en- 
ergy parameters by 10% or more. 

As shown in Figure [5j to recover an unbiased best 
fit value for dark energy parameters (in our case the 
fiducial model, wq = — 1 and wi = 0) the correct re- 
combination model must be used and hence one has 
to introduce an additional parameter, e a in the anal- 
ysis. Introducing this extra degree of freedom could 
in principle lead to an increase in the uncertainties 
on others parameters. We have verified that the in- 
creased uncertainty in the dark energy parameters is 
small. Including e a in the analysis increase lcr errors 
on these parameters less than a 10% percent yield- 
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ing a Wo — 0.12 and o~ wl = 0.15 (compare with Hill) . 
Forthcoming data sets will thus disentangle effects of 
non-standard recombination from the effects of dark 
energy. 

V. CONCLUSIONS 

In this paper we have investigated how a biased de- 
termination of the sound horizon due to a incorrect 
recombination model affects cosmological parameters 
measurements from future BAO data. We have shown 
that assuming a standard recombination model when 
the true model is different can bias the measure of 
the CMB sound horizon and the best fit value of cos- 
mological parameters such as the Hubble parameter. 
This shift propagates in a shift on the values of Da(z) 
and H(z) and hence on the best fit values of wo and 
W\ determined from BAO surveys. We have shown 
that a deviation from standard recombination which 
is still allowed by current CMB data, propagates in a 
bias in wq and w\ comparable to their la statistical 
error from forthcoming surveys. An incorrect calibra- 
tion of sound horizon due to a wrong assumption for 
the recombination model could bias estimation of dark 
energy parameters by 10% or more. In this case future 
surveys may thus incorrectly suggest the presence of 
a redshift dependent dark energy component. 

To recover an unbiased best fit value for dark energy 
parameters the correct recombination model must be 



used or must be described by the set of cosmologi- 
cal parameters used in the analysis. We have con- 
centrated on the case of delayed recombination where 
extra sources of resonance (Lyman-a) photons mod- 
ify the evolution of the ionization fraction. This is 
well modeled by the addition of a single extra param- 
eter. We have employed a Fisher matrix formalism to 
forecast errors and parameters biases for forthcoming 
CMB (Planck) and BAO (BOSS, WFMOS) surveys. 
For this data set combination, introducing the addi- 
tional degree of freedom describing deviations from 
standard recombination does not significantly degrade 
the error-bars on dark energy parameters and yields 
unbiased estimates. 

This is due to the CMB-BAO complementarity: 
forthcoming data sets will thus disentangle effects of 
non-standard recombination from the effects of dark 
energy. 
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FIG. 4: Shifts on best fit of Da(z) and H(z) (solid line) for the mean redshift of each bin caused by a wrong determination 
of sound horizon when fixing e Q = in the hypothesis that true value is e a — 0.05. Dashed contours are the 1,2 a 
constraints obtained assuming the correct recombination model while solid contours assume standard recombination, 
The star is the correct fiducial model for Da (z) and H(z). 
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FIG. 5: Constraints at 68% c.l. on wo and Wi, assuming 
an incorrect recombination model (solid line) and with the 
correct model (dashed line). The true model is indicated 
with the star. 




